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We study integral representation of so-called c?- dimensional Catalan numbers Cd{n), 



defined by 
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.p=0 



{n+py. 



{dn)\, (i = 2, 3, . . ., n = 0, 1, . . .. We prove that the Cd{nys are 

the nth Hausdorff power moments of positive functions Wd{x) defined on x G [0, d'^]. 
We construct exact and explicit forms of Wd{x) and demonstrate that they can 
be expressed as combinations of rf — 1 hypergeometric functions of type d-iFd-2 
of argument xjd'^. These solutions are unique. We analyse them analytically and 
graphically A combinatorially relevant, specific extension of Cd{n) for d even in the 

d/2-l 

is analyzed along the same lines. 
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Multidimensional Catalan and related numbers as Hausdorff moments 

I. INTRODUCTION 

Amongst many existing generalizations of classical Catalan numbers C{n) = ^^ ( ^) , 
those that include the parameter that in certain sense can be associated with the spacial 
dimension d, are particularly interesting. They permit to extend to higher dimensions d > 2 
the notions of objects enumerated by C{n) in d = 2. We shall be concerned in this note with 
one of such generalizations, called (i-dimensional Catalan numberai"-, which are defined as: 



C,(n) 



'd-l , 

pi 



n 



.^=0^^^ + ^)' 



{dny., n = 0,l,...; d = 2,3,..., (1.1) 



which for d = 2 clearly reduce the conventional Catalan numbers C{n). The form of Eq. (11. ip 
guarantees that Crf(O) = 1 for all d. We shall refer to Sloane's Online Encyclopedia of 
Integer Sequences (OEIS)^ and quote initial terms, n = 0, 1, . . . , 7, of several sequences 
Cd{n), d = 2,3,4 and 5, along with the labelling of their entries in the OEIS: 

• ior d = 2: 1, 1, 2, 5, 14, 42, 132, . . . (A00108), which are the Catalan numbers, 

• for rf = 3: 1, 1, 5, 42, 462, 6006, 87516, . . . (A005789, A151334), 

• ioT d = 4: 1, 1, 14, 462, 24024, 1662804, 140229804, . . . (A005790), 

• for rf = 5: 1, 1, 42, 6006, 1662804, 701149020, 396499770810, . . . (A005791), 

• for general d, see A060854. 

The explicit form of Cd{n)^s permits one to immediately write down some of their charac- 
teristics. If A(/c, a) = |, ^ii, . . . , "^^"-^ denotes a special list of k elements, then the ordinary 
generating function (ogf) of Cd{n)^s can be written as: 

g{d, z) = Y^Cdin) z- = dFd-i (^2^^"^' ^^ J^'^) • (l'^) 

Similarly, the exponential generating function (egf) of Cd{n)^s reads: 

G{d,z) = Y,Cd{n) -,=.Fd(^ f ^' ^\ d'^z) . (1.3) 

The use of Stirling's formula gives the leading term of n — )■ oo asymptotics for Cd(n): 

Cd(n) ,,^ n-'^rf'^" + . . . , rf = 2,3,.... (1.4) 

In Eqs. (II. 2p and (II. 3p we have employed the standard notation for the generalized hyper- 
geometric function pFg I ^Z'''^ x), with (a^) and (/3g) the lists of p "upper" and q "lower" 
parameters, respectively, see^. Observe that since in d^d of Eq. (II. 3p there is a pair of lower 
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and upper parameters differing by one, the appropriate function ^Fa can be reduced to a 
combination of d-i^d-is, see the formula 7.2.13.17 on page 439 o^. 

Inspired by the very fruitful interpretation of Catalan numbers C{n) as moments of a 
positive function on x G [0, 4], which is intimately related to the famous Wigner's semicircle 
law^, we set out to consider the sequences Cd{n), rf > 2 as Hausdorff power moments and 
have defined an objective of obtaining for d > 2 the equivalents of the solution for d = 2, 
quoted in Eq. (12.81) below. 

The paper is organised as follows: in Sec. Ullwe describe the method of obtaining exact 
and explicit solutions for d > 2. Subsequently we write down the general solution for d 
arbitrary and quote the specific cases of c? = 2,3,4 and 5. In Section UTTl we discuss some 
possible generalizations that has not been yet carried out. In Sec. [IV]we close the note with 
short conclusions. 

II. SOLUTIONS OF THE HAUSDORFF MOMENT PROBLEM 

A. We are seeking the solutions of the following Hausdorff moment problem: 

"Rid) 

x''Wd{x)dx = Cd{n), n = 0,l,..., d = 2,3,..., (2.1) 

where R{d) - the upper edge of the support of Wd{x) - will be determined below. The 
conventional estimate R{d) = lim [Cd(n)]^/" = d^ will be confirmed later by the Mellin 

n—>oo 

transform analysis. As a preliminary step we shall demonstrate that the sought for Wd{x) 
defined in Eq. (12. ip is positive. By using the Gauss- Legendre multiplication formula for 
gamma function to Eq. (II. ip and introducing complex s such that n = s — 1 we obtain 



/d-i \ d-ir s-l + i±i 

Cd{s - 1) = (2vr)^rf^^ mkl] {d'^) ^ n \. ... ^ (2.2) 

\fc=o / j=o y^^Ji 

which should be interpreted as the Mellin transform of Wd{x) i. e. j^ x^~^Wd{x)dx, denoted 
by M. ^d{x)'i s], see Ref.-. Since j > "'"'"^~ for all < j < rf — 1, the individual term labelled 
by j in the second product of Eq. (12. 2p has the inverse Mellin transfornt^.^ see the formula 
8.4.2.3 on page 631 front^ 



M 



-1 



r(.-i + ^ 



n-s+j) 



X 



^{j+l)/d~l Q _ ^y~(j+i)/d 

7 T— \ ) (2-3) 



j = 0,1, . . .d — 1, e. g. it is proportional to the standard probabilistic beta distribution^ in 
the variable x, which is a positive and absolutely continuous function for < x < 1. We 
perceive now Cd{s — 1) as a product of d such individual terms. Then the weight Wd{x) is a 
positive and absolutely continuous function on [0, R{d)], since it is a ci-fold Mellin convolution 
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of positive and absolutely continuous functions on [0, 1]. In the final result we accommodate 
the prefactor (t/'^)* which indicates, via elementary property of the Mellin transfornt^, that 
the solution of Eq. (12.11) will depend on xjd'^. 

It turns out that such a d-fold Mellin convolution can be carried out explicitly. The key 
step is first to identify the weight Wd[x) as a special case of Meijer G function G!^'"-. This 
is a direct consequence of (12. 2 p and reads 



Wa{x) = (27r)^ 



'^*-(n«)«2;S(|. 



0,l,...,ci-l 
-A(rf,0) 



(2.4) 



where A(n, a) = -, ^^, . . . , s+Hid, Next, the Meijer G function is converted to the hyperge- 
ometric form, using formulas 16.17.2 and 17.17.3 of^, which is the Slater theorem. We quote 
only the final result which is equal to: 



i=i 



I —1 — I 

d' d' ■ ■ ■ 

i-1 

d' ■ ■ ■ ' -^ d ' -^ ' d' 



^ + 2-i 



l_ll_2 l_i^.l,il,2 ^ 



d-j-1 
d 



X 
d^ 



(2.5) 

defined ioi < x < d'^ , which implies R{d) = d'^ in Eq. (12. ip . (For the reader's convenience 
we point out that in Eq. (12. 5p . in the lower list of parameters of d-iFd-2, there are two 
sequences of numbers, which contain j — 1 and d — l—j terms, respectively). The numerical 
coefficient Cj{d) is equal to 



Cj{d) = (27r) VrfJ-'^+l 



pi 



iiTip+j 



nr(i 

.fc=l 



d-1 



nr^ 



U=i+i 



(2.6) 



where j = 1, . . . ,d — 1 and d = 2,3, . . . . 

The structure of parameter list of Meijer G function in Eq. (12. 4p warrants that the 
assumptions of formula 2.24.2.1 in^ are satisfied: 

d-l d-1 



sE^-E^- 



d' 



<0, d = 2,3, 



(2.7) 



A:=0 



fc=0 



Therefore the Mellin transform of Wd{x) is well defined for ^{s) > ^-p. 



We shall explicitly write down the solutions for c? = 2, 3, 4 and 5, starting with W2{x): 

1 IA- X 



Woix) 



< z ^ 4, 



(2. 



27r V X 

which is obtained in many references^i^^, see Fig. [1] It is the only density that can be 
expressed by an elementary function. Furthermore, for d > 2 no density can be expressed 
by standard special functions, and the hypergeometric form is the final one. For li = 3, . . . , 5 
the solutions read: 



W.ix) 



ci(3) 



X 



1/3 



.Fl 



4 _1 
3' 3 

4 

3 



33 ] + ^2/3 2^1 



3' 3 
2 
3 



X 



— 1 , < a; ^ 33 



(2.9) 
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WJx) 



ci(4) 



X 



1/4 



3-^2 



C3(4) 



X' 



3/4 



3-^2 



9 5 


1 


4' 4' 


4 


5 3 

4' 2 




11 7 


3 



X 



4'' + 



4 ' 4' 4 
1 3 

2' 4 



X 



C2(4) / 


5 3 

2' 2' 


1 


3 5 

4' 4 





^h o<x^4^ 



(2.10) 



P^.fa;) 



ci(5) 



X 



1/5 



4-^3 



+ 



C3(5) 



X 



3/5 



4-^3 



< X ^ 5" 



16 _11 _ 

' 5 ' 5 ' 

6 7 8 
5' 5' 5 



5' 



18 _13 _8 
5 ' 5 ' 5' 
3 4 6 
5' 5' 5 



X 

X 
55 



C2(5) 



X 



2/5 



4-f'3 



, C4(5) 



X 



4/5 



17 


12 


7 


2 




5 ' 


5 ' 


5' 


5 






4 6 


7 




55 




5' 5 


5 




19 


14 


9 


4 




5 ' 


5 ' 


5' 


5 


X 




2 3 


4 




55 




5' 5' 


b 







(2.11) 



The coefficients Cj{d), j = 1, . . . ,d — 1, ioi d = 3, . . . ,6, are collected in Tab. [H With Cj(6)'s 
given there and using Eqs. (12. 5p and (12.61) . the reader can easily reconstruct Wq{x), which 
will not be reproduced here. The solution Ws^x) is represented in Fig. [2l 



c,(3) 

c,(4) 

c,(6) 



4^2^/3 



— r, - 

10^ \ 3 



757r4 



43 



5Sv^r(|)'^4 59^r(f)V 



253^11^7r5B 



2"3^'' E 
7"'133l8057r"^ 



267317^5^1 



3i«D__ 
2i"73657rB 



4851-^ \, 4 



5''V5(AB)-i 59V5(AB)-^ 



2'^34l32r 



22\ 



2'^34 72l9r 



3I" 
11^53 2«7D 



2^^ 3^; 
ll45-*17-*23^29£; 



TABLE I. The coefficients Cj{d), see Eqs. (|2.5p and (|2.6p for d = 3, . . . , 6 and j = 1, . . . ,d — 1; 
To simplify the notation in Cj{5) and Cj{6) we set A = sin ( f ) , -B = sin ( ^ j and £■ = r( | j , 



III. GENERALIZATION OF MULTIDIMENSIONAL CATALAN 
NUMBERS 



In this paragraph we analyse the extension of C^^n) obtained by replacing (4n)! in 



.r=0 



Eq. (II. ip by (2n)!(2n+2)! . The corresponding sequence Di{n) = 6(2n)!(2n+2)! 

has attracted attention in several contexts, as it appears in^*^"— . 

The initial terms of D^^n) are 1,1,4, 30, 330, 4719, 81796, 1643356, forn = 0, 1, . . . , 7. It 
is listed as A006149 in OEIS where also additional information can be found. It turns out 
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W^(x) 0.4 




FIG. 1. The density W2{x), s. Eq. (fTg]) . 
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FIG. 2. The density W'iix), s. Eq. ^^. 



that the ogf of D^iji) can be expressed by the elliptic functions lE(y) and ]K(y)^: 



n=0 



1 + 62 n - 162) d + 112^) 



An)z- 



4^2 



240 TT z^ 



K(4v^) 



1 + 224^ + 256^2 

240 TT ^3 



E(4v^). (3.1) 



In fact the sequence D^iji) allows for the same kind of analysis as does the ensemble of 
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Cd(n)'s. The Hausdorff moment problem for D^{n), namely 

\^ V,{x) dx = D,{n) = 6(y (^^ + ^)' ^ n = 0,l,..., (3.2) 

r=0 

can be solved by the method of Mellin convolution and the use of Meijer G function elu- 
cidated above. The weight can be proven to be positive on x G [0, h] with h = 16 and 
reads: 



VJx) 



157r2 



64 



+ 5Q^/x ■ 



,.3/2 



EL/1 



16 



2v^(16 + x)kL/i 



16 



(3.3) 



We plot the function V4(x) on Fig. |3l 



VAX) 




10 12 14 16 



FIG. 3. The density Vi{x), s. Eq. iK3\i . 



The sequence D^{n) analysed above is a special case rf = 4 of the following generalization 
of Cd(n) defined for even d: 



DAn) 



'd-i 

n 

r=0 



{n + r)! 



d/2-l 

n 

s=0 



(2n + 2s)! 
(2^)! 



ci = 2,4,6 



, i, v^, . . . , 



(3.4) 



in which the parameter d should not be associated anymore with the spacial dimension. 
Several exact characteristics of the sequences Dci{n) are available. The ordinary generating 
function reads 

1 3 d-l 1 

(3.5) 



^(^, . )=./.i^/.-i (,--,- 2^. :.,rf 



2'^z 
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0.0020 



0.0015 



Ve^x^ 0.0010 



0.0005 




FIG. 4. The density Veix), s. Eq. ([SSD for d = 6. 



whereas the corresponding exponential generating function is equal to 



G{d,z)=,F, 



1 3 



d-1 



1,2 



2! 2' ■ ■ ■ ' 2 ' -^' "'•••' 2 

1,2,3,. ..,d 



2'^z 



(3.6) 



The leading term of the n — )■ oo asymptotics for Dd{n) can be obtained by using the Stirling's 
formula and it has the following form 



It is remarkable that the Hausdorff moment problem for Dct{n), i. e. 

r-K(d) 



(3.7) 



x^VAx)dx = DAn), ?T. = 0, 1 



, J-, . . . , 



d = 4,6. 



can be exactly solved as well in terms of positive functions Vd{x) defined on x G [0, 2°'] 
i. e. K{d) = 2'^, which read: 



d-l 



VJx) 



l,...,d-l 



r=0 ^d/2,0 I j^i 2' 2 ^ 

d/2-1 , . '^d/2,d/2 I Ofi' _1 1 3 d^ 

n T(k+'.)k\ V 2'2'2'---' 2 



(3.^ 



fc=0 



Here the condition 2.24.2.1 in^ implies — |((i+ 1) < 0, d = 4, 6, . . ., which is always satisfied. 
In addition, the Mellin transform of Vd{x) is well defined for 9fJ(s) > |. The proof of positivity 
of Vd{x) can be carried out along the lines exposed in the previous Section. 
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Since in Eqs. (13. 8 p in both parameter lists in the Meijer G function there are index 
pairs that differ by an integer, this Meijer G function cannot be represented by a sum of 
generalized hypergeometric functions. However the expression (13. 8 p can be easily manipu- 
lated algebraically and represented graphically—. In Fig. |l]we display Vq{x) in the range 
X G [10, 60]. Observe the rapid decrease of this function for x > 25, followed by a large 
region where it is practically flat and equal to zero. Similar behavior is observed for higher 
values of d. 

IV. DISCUSSION AND CONCLUSIONS 

We have treated in this work essentially two generalizations of conventional Catalan 
numbers, which are related to such notions as Young tableaux, hook lengths, generalized 
Dyck paths, etc.—. They all turn out to be moments of positive functions on supports 
included in the positive half line. The relevant weight functions have been obtained explicitly 
and analyzed graphically. All these positive functions are unique solutions of Hausdorff 
moment problems. The key tool in this approach had been the inverse Mellin transform 
and the encoding with Meijer G functions. The positivity of solutions has been rigorously 
proven using the method of Mellin convolution, applied to related problems previously^"—. 
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